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Abstract 

Gaussian linking of superconducting loops containing Josephson junctions with enclosed mag- 
netic fields give rise to interference shifts in the phase that modulates the current carried through 
the loop, proportional to the magnitude of the enclosed flux. We generalize these results to higher 
order linking of a superconducting loop with several magnetic solenoids, and show there may be 
interference shifts proportional to the product of two or more fluxes. 
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I. INTRODUCTION 



Interference effects, both constructive and destructive, are mainstays for distinguishing 
between quantum and classical phenomena. Examples include interfering scattering am- 
plitudes, both bosonic and fermionic, formation of condensates, entanglement, etc. The 
Aharonov-Bohm effect describes the self-interference of a charged particle that can 
travel along two semiclassical paths whose combined path is gaussian linked with a mag- 
netic solenoid carrying flux $. The measurable phase shift is oc $. We have argued in 
Ref. that there could exist generalizations to cases of higher order linkings. The simplest 
example is a Borromean ring arrangement where the semiclassical path corresponds to one 
ring, which has higher order linking with two flux tubes carrying fluxes $i and $2? which 
make up the other two rings. We found the phase shift in this system is oc $i$2- Higher 
order cases were explored in Ref. ^, 4| and shown to be related to commutator algebras of 
homotopy generators of the configuration space M'^\{Ti UT2}, where Ti and T2 are the tubes 
containing the fluxes. The same general logic can be applied to systems of superconductors, 
Josephson junctions, and magnetic fluxes where the Josephson effect can arise Here we 
will study interference in a symmetric arrangement of two identical semicircular supercon- 
ductors joined by two identical Josephson junctions and derive the response of such systems. 
We conclude with a discussion of possible applications. 

In the case of gaussian linking of a loop of superconductor with a magnetic solenoid, the 
Mercereau effect |6|] is due to the phase change in the macroscopic wave function, which is 
in turn related to the currents in the superconducting components. The effect is due to the 
presence of a vector potential A, which is the fundamental object responsible for the phase 
change. Exploration of higher order linking is again due to the presence of a vector potential 
but in these instances it requires careful choices of gauge. 



II. THE JOSEPHSON EFFECT 



It will be sufficient for our purposes to consider a macroscopic model of superconduc- 
tors. Following Feynman we approximate the superconductors coupled via a Josephson 
junction as a two-level system. Let tpi and tp2 be the states, Ei and E2 the energy lev- 
els of the superconductors. The Schrodinger equation for the coupled system of the two 
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superconductors becomes 



ih{dtPi/dt) = Ei^i + K exp {i^)^2, (1) 
ih{d^2/dt) = E2^p2 + exp (2) 

where K is the couphng energy and (pis a. phase, which arises from the most general hermitian 
hamiltonian 2x2 matrix. The dependence of (p on the vector potential A can be found 
from gauge invariance considerations. For a gauge transformation A A + V/, ipi i— >■ 
ipi exp {iqf /2hc), 1^2 ^ exp {—iqf/2hc) with an arbitrary function / of space coordinates, 
we find (p ^ (p + qf /he, from which it follows that (p = {q/hc) J A ■ dx. Here g = 2e is the 
charge of an electron pair. 

After the substitutions 'ipi = iV'il exp {i9i) and V'2 = \4'2\ exp (2^2), the Schrodinger equa- 
tion becomes 

h{d\iPi\ydt) =2K\tP^\\i;2\sme, (3) 
h{d\ij2\Vdt) = -2K\tlji\\tlj2\sme, (4) 
-h\iJi\{d9i/dt) = EilV^il + K\ij2\ cos9, (5) 
-h\Hide2/dt) = E2\iJ2\ + K\tP^\ cos9, (6) 

where 6 = (p + 62 — Oi- The current from superconductor 1 to superconductor 2, 
which is equal to minus the current from superconductor 2 to superconductor 1, is thus 
/ = {2K/h)\ilJi\\tp2\ sinO. (In a self-consistent computation, a current from a battery which 
connects the two superconductors is also included. The result for the superconducting 
current is precisely J; see, for example, Ref. j^.) The electron densities in the two supercon- 
ductors are approximately equal and independent of time; let p be this common constant. 
This gives I = Josin^, where Iq = 2Kp/h. Integrating the phase equations, we find 

e{t) =(P + 6(0) + h-^ [ dt' (Ei(t') - E2{t')). (7) 

Jo 

The quantity (-Ei(t) — -E2(i!^))/Q_ 'represents an electric potential applied to the junction. The 
dc and ac Josephson effects [5| arise for \Ei{t) — E2{t)\ <^ K and \Ei{t) — -£^2(^)1 ^ K, 
respectively. 

Our interest here is in the Josephson effect with zero potential across the junction, Ei{t) — 
E2{t) = 0, and nonzero magnetic field constrained to the opening of the superconducting 
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FIG. 1: A diagram of an experimental setup for the detection of the Josephson effect. C and 
C" are paths from the point P to the point Q through the superconductors with the Josephson 
junctions J' and J" and the total current / from P io Q. Ci is the magnetic solenoid carrying flux 
$1. The Josephson effect (for a review see gaussian) linking of 

the closed curves C = C'C"~^ and Ci. 

ring with two Josephson junctions; see Fig. [H Let 6' and 6" be the phase changes due to 
the vector potential Ai of the currents through the junctions J' and J". The phase changes 
from the point P to the point Q along the paths C and C" are 

cp' = e' + (q/hc) [ Ai-dx, (8) 
Jc 

(J)" = e" + {q/hc) [ Ai-dx. (9) 
JC" 

Since the wave function is single valued, this requires 0' = 0", and so we find 9" — 9' = 
27r$i/<l'o- Here $i = §^ Ai ■ dx is the fiux due to the solenoid along Ci passing through a 
surface spanned by a closed curve C = C'C"^^ and $o = 27r^c/g is the fiux quantum. The 
total current from the point P to the point Q is 

/ = /o sin {\{9' + 9")) cos (7r$i/<l>o). (10) 

For a fixed value of $i, the corresponding maximal total current is 

/max = /o|cOs(7r<l>i/$o)|, (H) 

which itself has maxima when $1 = n$0; n eTj. 

The fiux is actually $1 = $i,cxt + LI where $i^ext is the external fiux through the loop, 
L is the self-inductance, but here and in what follows we assume L is negligible. (We have 



made a number of simplifying assumptions, for example, that self-inductance of SQUID 
components are negligible, none of which, if relaxed, affect our basic conclusions.) 

We will call the phenomena reviewed in this section first order Josephson effects to dis- 
tinguish them from their generalizations which we now proceed to describe. 

III. THE SECOND AND HIGHER ORDER JOSEPHSON EFFECTS 




FIG. 2: A diagram of an experimental setup for the detection of the second order Josephson effect. 
Ci and C2 are the magnetic solenoids. C and C" are paths from the point P to the point Q 
through the superconductors connected by Josephson junctions J' and J" . The total current from 
P to Q is /. Ci and C2 are magnetic solenoids carrying fluxes $1 and <I>2- The second order 
Josephson effect is due to the second order linking of the set of three closed curves C = C'C"~^, 
Ci and C2. 

Now consider the case where we have two solenoids carrying magnetic fluxes $1 and $2 
and whose center hues run along Ci and C2, and a superconducting ring along the closed 
curve C = C'C"^^ with two Josephson junctions J' and J" in parallel as shown in Fig.[2l The 
two solenoids and the superconducting ring are in a Borromean rings configuration. Note 
that in this arrangement neither Ci nor C2 has gaussian linking with the superconducting 
ring C, nor do Ci and C2 link with each other. However, the set of three rings C, Ci,C2 
is indeed linked. This second order linking and its higher order generalizations is what 
will lead to our results. In other words, we will find that even though our system lacks 
first order (gaussian) linking, a phase difference can still exist upon traveling around the 
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superconductor. To find this phase, we must choose a gauge that detects it. Such a gauge is 
[2! A12 = 1^2(71^2 — 72^1)- Here Ai and A2 are the vector potentials due to the solenoids 
along Ci and C2, the quantities 71 and 72 are defined by 

7i = Sj + (q/hc) j Aj ■ dx, (12) 

where 61 and 62 are constants, and F is a path that runs along C. The quantity k2 is a 
normalization constant, the value of which we discuss below. 

Using A12 in the phase integral and following computations in Ref. j^, we find 6" — 6' = 
47r^fc2$i$2/$o- The total current from point P to point Q is 

/ = Jo sin {\{d' + 6")) cos (27r2A;2$i$2/$^). (13) 

For fixed values of $1 and $2, the maximal total current flowing in the superconductor is 

/max = h I COS {2t:'' k2^ 2 1 ^l) | ■ (14) 

The smallest value of the constant A;2 > for which the fluxes $i = mi$0) "^2 = '^2*&0; 
where mi, 1712 G Z, lead to maxima of the quantity /max is k2 = (27r)^^. This is precisely 
the value we obtained in Ref. ^ by imposing an analog of the Dirac string condition on 
the second order phase for the Aharonov-Bohm effect. Nevertheless, the value of ^2 must 
ultimately be determined by experiment. 

Also, for the value k2 = (27r)^^, if either $1 or $2 is equal to $0 or — $0) then in terms of 
the other flux, appropriately relabeled, the expression (fT4|) for the second order /max reduces 
to the expression ( fTTl) for the flrst order /max- 

An essential feature of the above result is that the current / is a periodic function of 
the quantity 7r/c2$i$2/^i..^with period 1. We can also derive this property by modifying the 
method used by Block ll| for the flrst order Josephson effect as follows. 

The total gauge potential includes internal and external parts, A = Am + Agxt, the 
external magnetic fleld being due to the external sources. Assuming that the external fleld 
V X Aext vanishes inside the superconductors, we can write Aext = V7ext- As a result, A 
is a gauge transformation of Am, and so 

?/'(Am + Acxt) = ^(Ain) exp (ig7ext/M- (15) 

Since ip{A) is single valued, we flnd that ipi-^m) is multiplied by the factor 
exp (— 4ot^A;2$i$2/^o) ^-fter the charge q travels around a closed curve C. This factor is a 



periodic function of 7rfc2 $1^2/^0 with period 2. This imphes the same periodicity property 
for the wave function ?/'(Ain) and the energy E. Assuming time reversal symmetry as in 
Ref. we find that the free energy and thus the current, which is given by minus the 
derivative of the free energy with respect to the external flux, are both periodic functions of 
7rA;2$i$2/*^'o with period 1, in agreement with the result proved earlier. 

More generally, it is straightforward to arrange n solenoids with the fluxes $1, . . . , 
and a superconducting ring in such a way that they are linked with nonzero nth order 



linking 10|]. We similarly find the phase difference 

r-^' = (27r)"A;„$i---$„/<l>^, (16) 

the current 

/ = Jo sin iliO' + 6")) cos (|(27r)"A;„<l>i ■ ■ ■ $n/$^), (17) 
and its maximal value for fixed values of $1, . . . , $„, 

/max = /ol COS (i (27r)" A;„<l>i • ■ • $n/$o) I- (18) 

Other properties of these systems can be investigated. 

The smallest value of the constant A;„ > for which the fluxes $j = rrij^o, where rrij G Z, 
lead to maxima of the quantity /max is kn = (27r)^^". Again this is precisely the value we 
obtained in Refs. 3, ^ by imposing an analog of the Dirac string condition on the phase for 
the nth order Aharonov-Bohm effect. Nevertheless, as pointed out with the k2 case above, 
the value of kn must ultimately be determined by experiment. 

Also, for the value /c„ = (27r)^~", if one of the fluxes is equal to $o or —^q, then in terms 
of the remaining fluxes, appropriately relabeled, the expression f|T8l) for the nth order /max 
reduces to the analogous expression ( JTSl) for the (n — l)st order /max- 

Similar to the case n = 2 above, we can modifying the method used by Block and prove 
that for any n the current / is a periodic function of the quantity |(27r)"^^A;„$i ■ ■ ■ $„/$q 
with period 1. 

IV. CONCLUSION 

We generalize the Josephson effect to its higher order analogs in which a superconducting 
loop links with several magnetic solenoids and the resulting interference shifts are propor- 
tional to the product of two or more fluxes. 



One can conceive of a number of applications for devices build to take advantage of higher 
order linking. Such a system could be less invasive than first order devices because it could 
keep the SQUID some distance from an experimental sample. Possible applications include 
both rf and dc SQUlDs that measure higher order linking of multiple fluxes. Under some 
circumstances such devices could be useful in measurements of complex biological systems, 
or any systems where direct gaussian linking of a magnetic flux with a SQUID is impractical, 
but where higher order linking is possible. For example, one could have a system of (i) a 
fixed but adjustable flux tube, i.e., a solenoid; (ii) an unknown flux to be measured, and (iii) 
a SQUID. If the three components can be arranged to have higher-order linking, then the 
unknown flux could be measured, even though it has no gaussian linking with the SQUID. 
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